The VLT-FLAMES survey of massive stars: rotation and nitrogen 
enrichment as the key to understanding massive star evolution 

1. Hunteri'2, I. Brott^ D.J. Lennon^ N. Langer^, P.L. Dufton\ C. Trundle^ S.J. Smartt^ 

A. de Koter^ C.J. Evans^ and R.S.I. Ryans^ 

ABSTRACT 

Rotation has become an important element in evolutionary models of mas- 
sive stars, specifically via the prediction of rotational mixing. Here, we study 
a sample of stars, including rapid rotators, to constrain such models and use 
nitrogen enrichments as a probe of the mixing process. Chemical compositions 
(C, N, O, Mg and Si) have been estimated for 135 early B-type stars in the Large 
Magellanic Cloud with projected rotational velocities up to ~300kms~^ using a 
non-LTE TLUSTY model atmosphere grid. Evolutionary models, including rota- 
tional mixing, have been generated attempting to reproduce these observations 
by adjusting the overshooting and rotational mixing parameters and produce 
reasonable agreement with 60% of our core hydrogen burning sample. We find 
(excluding known binaries) a significant population of highly nitrogen enriched 
intrinsic slow rotators [vsini <50kms~^) incompatible with our models (~20% 
of the sample). Furthermore, while we find fast rotators with enrichments in 
agreement with the models, the observation of evolved [logg < 3.7dex) fast ro- 
tators that are relatively unenriched (a further ~20% of the sample) challenges 
the concept of rotational mixing. We also find that 70% of our blue supergiant 
sample cannot have evolved directly from the hydrogen burning main-sequence. 
We are left with a picture where invoking binarity and perhaps fossil magnetic 
fields are required to understand the surface properties of a population of massive 
main sequence stars. 
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Introduction 



Stellar rotation in massive stars has been adopted in recent theoretica " 



ample, predict the correct blue to red supergiant ratio flMaeder fc Mevnet 



itors of gamma-ray bursts through homogeneous evolution flYoon Sz Langerll2005[) and Wolf- 



models to, for ex- 



2005 



20011). t he progen- 



Vink fc de Koter 



Raye t populations as a function of metallicity (iMeynet fc Maedei 
20051 ). However, the consequent surface enrichment of helium and nitro gen through rota- 
tional mixing is poorly constrained by observations (jPaflon et al.ll200ll ). with the mixing 
typically c alibrated from eit her studies of main-sequence stars with low projected rotational 
velocities (IKorn et al.ll2002l ) or evolved supergiant stars (jVennlll999l ). The former are biased 
towards slow rotators while the latter have evolved beyond the core hydrogen burning phase. 
The analysis of an unbiased sample of fast rotating core-hydrogen burning stars in order to 
properly calibrate the predicted rotational mixing efficiency was a strong drive r for a large 



survey of O - and early B-type star s in o ur Galaxy and the M agellanic Clouds (lEvans et al. 



2005. 200fil ). iHunter et all (l2007al ) and iTrundle et al.l (120071 ). hereafter Paper I and II re- 
spectively, have presented chemical abundances for the narrow lined objects (predominately 
slow rotators) in the LMC sample. Here we extend the study with the chemical analysis of 
the faster rotating stars. 



2. Observations and model atmosphere analysis 



High-resolution [R ~ 20 000) spectra from the Fibre Large Array Multi-Element Spec- 
trograph (FLAMES) at the European Southern Observatory Very Large Telescope were 
obtained for some 750 O- and early B-type stars located towards clusters in our Galaxy 
and the Large and Small Magellanic Clouds (LMC and SMC respectively). A discussion of 



target select ion, observational details and initial data reduction can be found in lEvans et al. 



( 12005. 20061 ). In this letter we discuss the LMC early B-type stars, using nitrogen as a probe 
of chemical enrichment. 



A grid of non-LTE TLUSTY model atmospheres (IHubeny fc Lanzlll995l ) has been used 
to calculate the atmospheric parameters and chemical cornposition s (C, N, O, Mg and Si) of 
our targets as described in Papers I, II and iHunter et al. J2007bh . hereafter Paper HI. We 
have fitted the lines with rotationally broadened profiles to estimate the equivalent widths 



3 



(EWs), since at significant rotational velocities (>50kms~^) the line shape is rotationally 
dominated. For those objects in which no nitrogen features could be identified, an upper 
limit to the nitrogen abundance was estimated by placing an upper limit on the equivalent 
width of the strongest Nil line in our spectral range, which is located at 3995A (Paper I). 

The mean abundances of C, O, Mg and Si are in excellent agreement (within 0.05 dex) 
with the LMC baseline abundances given in Paper I. However the mean nitrogen abundance is 
0.36 dex higher than its baseline abundance indicating that nitrogen enrichment has occurred 
in many of the stars. Additionally the scatter in the nitrogen abundances is over a factor 
of two larger than that of the other elements, indicating that different levels of enrichment 
explain the mean nitrogen abundance, rather than systematic errors. Since nitrogen is an 
important element in the CNO-cycle, the surface nitrogen enrichments can be used as a 
measure of the mixing efficiency. Note that although a corresponding carbon depletion of up 
to ~0.2dex would be expected, within the uncertainties in determining carbon abundances 
(see Paper I), such an effect would be difficult to observe. In Fig. [1] the nitrogen abundances 
are plotted as a function of the projected rotational velocity with the sample being split into 
two groups; core hydrogen burning objects with surface gravities >3.2dex and supergiants 
having surface gravities <3.2dex (see Paper III). 



3. Stellar evolution models 



Comparison with published stellar evolution models such as the Geneva models (IMaeder fc Meynet 



200ll ) is complicated as no rotating models at the LMC metallicity are currently available. 
Additionally, the initial chemical composition for evolutionary models is often scaled from 
solar composition, and this is not appropriate for all elements, in particular nitrogen. 

New evoluti onary model s have been calculated (Brott e t al. in prep.) usin g the stellar 
evolution code of lYoon et al.l (120061 ) . which i ncludes rot at ion (jHeger et al.ll2000l ) and angular 
momentum transport by magnetic torques (jSpruitll2002l ). Only two differences ap ply here: 
we updated the mass loss prescription and now use the recipe of Vink et al.l (120011) . and we 



2002), which — in 



disregard the magnetically induced chemical diffusion term of ( SpruitI 
contrast to the magnetic angular mo nientuni tran sport — is not observationally supported, 
and appears controvercial at present (ISpruitll2006l ). 



As initial composition we adopted the LMC C, N, O, Mg and Si abundances given in 
P aper I and all the othe r metal abundances decreased by 0.4 dex from the solar composition 



of I Asp fund et a. 



( iPeimbert et al. 



Jj 20051) . Based on the recent discussion of the primordial helium abundance 



20071 ). we have updated the initial helium mass fraction for our LMC models 
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to Y= 25.6%, which together with the metalhcity of the adopted chemical composition (Z= 
0.5%), results in a hydrogen mass fraction of X= 73.9%. 

In Paper III we determined that the end of core hydrogen burning occurs at a logarithmic 
surface gravity of ~3.2dex. We find that a convective core overshooting of 0.335 pressure 
scale heights is required to reproduce this result. While this value is larger than what 
is typically assumed, with consequences for the post-main sequence evolution which still 
remain to be investigated, it provides the only way to understand the sharp drop in the 
rotational velocity distribution of our sample at a surface gravity of 3.2 dex. Evolutionary 
models neglecting overshooting indicate that the surface nitrogen enrichment is smaller by 
less than 0.15 dex at core hydrogen exhaustion for a rapid rotator. Hence, our principal 
conclusions are not significantly affected by the adopted amount of overshooting. 

The efficiency of rotationally induced mixing in our models is controlled by the pa- 
rameter /c, which is the ratio of the turbulent viscosity to the diffusion coefficient that 
describe s the transport of angular momentum by rotationally induced hydrodynamic insta- 
bilities ( iHeger et al.ll2000l ). The mean projected rotational velocity of the stars shown in 
Fig. [T] is llOkms"^, with the mean surface nitrogen abundance being 7.2 dex. Although 
this mean p rojected rotational ve locity is lower than that observed for cluster stars (see, 
for example, iHuang fc Giesll2006al ) there is no bias towards stars with low rotational veloc- 
ities. Our stars are a relatively unbiased sample of the true populations, as explained in 
Evans et al. I (120061 ). This population is dominated by field stars, which are known to rotate 
slower than cluste r stars , and our velocity distribution is comparable to the Galactic field 
stars of lAbt et all (120021 ). Additionally Be-type stars are excluded, hence our mean velocity 
is representative only for normal B-type stars and there is no velocity bias for these B-type 
stars. 



The parameter has been calibrated to reproduce the mean surface nitrogen abundance 
of the core-hydrogen burning stars at core hydrogen exhaustion for a 13 M0 model (the mean 
mass of our non-supergiant stars) initially rotating at 140kms~^ (110x4/7r to account for 
random an gles of inclinatioii ) . We find /c=2.28xl0~^, compared to 3.33x10"^ adopted 
previously (IHeger et al.ll2000l ). Using these parameters, we computed a grid of models for 
masses representative of our sample and a range of initial rotational velocities. 



4. Discussion and conclusions 



In Fig. [T]the nitrogen abundances are plotted against the projected rotational velocities 
for our sample stars, and compared with the evolutionary models. Given the large number of 
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objects in our sample, it is reasonable to assume a random distribution of inclination angles, 
and hence the rotational velocity of the tracks has been scaled by 7r/4. Although we have 
attempted to contrain the evolutionary tracks to the observations, we are able to reproduce 
agreement for only 60% of the data (excluding known radial velocity variables). 

However, two groups of core hydrogen burning stars in Fig. [T](a) stand out as being 
in conflict with the evolutionary models. Group 1 contains fast rotators which have little 
chemical mixing. In particular, this group includes fast rotating single stars with surface 
gravities indicating that they are near the end of core hydrogen burning (excluding radial 
velocity variables this is ~20% of the core hydrogen burning sample.) The gravities used 
to discriminate these stars have not been corrected for the fact that we may observed them 
almost equator-on, i.e. their true polar gravities could be larger, meaning that they are less 
evolved than our derived gravities i mply. However, applying the surface gravity corrections 



described in iHuang k. GiesI (j2006bl ) would increase the gravities by up to ~0.3dex for our 
faster rotators. As the zero age main-sequence gravity of our 13 models is about 4.3 dex, 
and since most of the stellar expansion on the main sequence occurs towards the end of core 
hydrogen burning, this gravity uncertainty can not reconcile the situation. 

The second discrepant group of stars in Fig. la are the 15 slow rotators [v sin i <50 km s~^) 
that show significant nitrogen enrichment (Group 2). This group also forms ~20% of the 
non-binary core hydrogen burning sample. Although it could be argued that they may be 
fast rotators observed pole-on, statistically this is unlikely. In order to reproduce their mean 
nitrogen abundance requires a rotational velocity of ~200kms~^. Using simple geometry to 
calculate the solid angle restriction that we must impose, i.e. 47r(l — cosO), and assuming 
that we need 15 stars to have this velocity, we expect less than one star to appear with such 
a low rotational velocity due to random inclinations. The simple corollary of this is that if 
we have 15 stars with velocities ~200kms^^ populating Group 2 through random angles of 
inclination, over 300 similary enriched stars with vsini values between 150 and 200kms~^ 
would be expected. Such a population clearly does not exist. 

Evolved rapid rotators with low nitrogen enrichment (Group 1) may be produced by 
close binaries with an initial period small enough to ensure tidal locking and slow rotation but 
large enough for highly non-conservative mass transfer, i.e. with an amount o f mass accreted 



which is sufficient to spin-up the star but insufficient to enrich it significantly (jPetrovic et al. 



20051 ). However, for many of the evolved stars in Group 1, there is no indication of binarity. 
The single star nature for these objects would pose a serious challenge to the theory of 
rotational mixing. 



While binaries m ight be able t o pop ulate Group 2 flLanger 1 12008| ) , another explanation 



may be more likely. iMorel et al.l (120061 ) have analyzed ten slowly rotating Galactic early 
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/?-Cephei B-type stars and found a highly enriched group (four out of ten) of which three 
have detected mag netic fields. A l thoug h pulsations in these stars may cause an apparent 
surface enrichment, iBourge et al.l (120071 ) report that only slight nitrogen enrichments would 
be expected. As such, assuming that the effect of pulsations is negligible, there appears to be 
a correlation between magnetic fields and nitrogen enrichment. Indeed, IWolff et al.l (120071 ) 
have attributed the large number of slow rotators often seen in massive star populations 
to magnetic l o cking, of the star to the accretion disk during the star formation process and 
Alecian et al.l (120071 ) also show a correlation between slow rotation and magnetic fields for 
magnetic A- and B-type stars. Hence we postulate that the highly enriched slow rotators in 
the LMC are analogs of the enriched Galactic magnetic stars. This would imply that, inde- 
pendent of metallicity, a significant fraction of early B-type stars are intrinsically magnetic, 
analogous to the well known situation for lower mass star s. In this context, th e identification 
of three He-rich slowly rotating OVz stars in NGC 346 (IMokiem et al.l |2006| ) suggests that 
the phenomenon of intrinsic magnetic fields in massive stars may not be confined to the 
B-star regime. If these mag netic fields are of fossil origin, i.e. possessin g long-term stability 
( iBraithwaite fc Spruitll2004l ). in contrast to the fields produced by the (ISpruitll2002l ) mecha- 
nism, one may speculate that stars in Group 2, and their Galactic counterparts, might be the 
progenitors of magnetars, analogous to the suggestion of iFerrario fc Wickramasinghd (120061 ) 
that magnetic white dwarfs evolve from magnetic A, F and late B stars. While the fossil 
field hypothesis might well explain the slow rotation, the physical process which leads to the 
enrichment of nitrogen in the stars of Group 2 remains to be identified. 

The nitrogen abundances of the supergiants (Fig. [^b)) appear to fall into two dis- 
tinct groups, one group having a level of nitrogen enhancement consistent with that seen 
for the majority of the core hydrogen burning objects (<7.2dex; Group 3), and a second 
having a much greater level of enrichment (>7.6dex; Group 4). It should be noted that 
although the Group 4 objects appear to be well fitted by evolutionary tracks with initial 
rotational velocities of 200-300 km s"^, such a rotational velocity distribution for supergiants 
is inconsistent with that for core-hydrogen burning objects and rotation does not increase 
the blue-supergiant lifetime. Hence, the enrichments of the supergiants in Group 4 (70% of 
the non-binary supergiant sample) are incompatible with the theory of rotational mixing. 
Additionally it should be noted that the lowest gravity stars (<2.8dex) are all enriched 
in nitrogen, implying that some evolutionary process not accounted for in the models is 
responsible for the enrichment. 

The simplest way to interprete this would be to characterize Group 3 as pre-red super- 
giant objects, and Group 4 as post-red supergiant objects. The nitrogen abundances and 
rotati onal velocities in Grou p 4 are indeed consistent with the predictions for a blue loop 
stage (iHeger fc Langerll2000l ). However, we note that our models do not return to the high 
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effective temperatures of the stars in Group 4 and hence their position on the H-R diagram 
cannot be reproduced. As seve ral of the enriched ob jects show evidence of binarity, mass 



transfer may also be important flWellstein et al.ll200ll ). 



To conclude, our study can not provide unambigous evidence for rotational mixing 
acting in massive stars. The fast rotators in our sample can be interpreted in two ways. 
Firstly as rotationally mixed stars, if the relatively non-enriched stars (Group 1) can be 
explained by binary effects. However, if it were confirmed that these Group 1 stars are single 
stars, then the enriched rapid rotators may need to be understood as binary products — and 
rotational mixing would be inefficient, or much more complex than described by the present- 
day shellular models. Additionally, the population of intrinsically slow rotators with nitrogen 
enhancements of up to a factor of ~6 implies that studies of nitrogen enhancements with 
stellar samples which are restricted to low projected rotational velocities — which applies 
to most previous works — are not suited for constraining, or motivating, the adoption of 
rotational mixing in massive stars. Finally our supergiant sample implies that these stars 
cannot be considered as representative of the amount of mixing that occurs during the 
hydrogen burning main-sequence and hence should not be used to constrain this process. In 
summary, our study provides a challenge to rotational mixing, and provides evidence for two 
other processes affecting the rotation and the surface abundances of populations of massive 
stars, likely binarity and magnetic fields. 
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Fig. 1. — Nitrogen abundance (12 + log [N/H]) against the projected rotational velocity 
(vsini) for core hydrogen burning (a) and supergiant (b) objects. Open symbols: radial 
velocity variables; downward arrows: abundance upper limits; dotted line: LMC baseline 
nitrogen abundance. The mean uncertainty in the nitrogen abundance is 0.25 dex while that 
in vsini is 10% and these are illustrated. These errors are largely independant of rotational 
velocity since the systematic uncertainties are comparable to the measurement errors. The 
bulk of the core hydrogen burning objects occupy a region at low vsini and show little 
or modest nitrogen enrichment. The tracks are computed for an initial mass of ISM© (a) 
and 19 (b), corresponding to the average mass of our non-supergiant and supergiant 
stars, respectively, and their rotational velocity has been multiplied by tt/A. Although the 
plot contains stars with a range of masses, comparison of the tracks shown in (a) and (b) 
show that any mass effect is negligible compared to the abundance uncertainties. In (a) the 
surface gravity has been used as indicator of the evolutionary status and the objects (see 
legend) and tracks have been split into red and blue to indicate younger and older stars 
respectively. However, this is illustrative only, since the evolutionary status is obviously 
continuous and not discrete. Gray shading in panel (a) highlights two groups of stars which 
remain unexplained by the stellar evolution tracks. In panel (b), gray shading highlights the 
apparent division of the supergiants into two distinct groups. The surface gravity estimates 
of many of the objects in Group 3, and the two apparantly rapidly rotating unenriched 
supergiants, are consistent with being in the core hydrogen phase within their uncertainties. 



